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p^' Abstract 

^^ \ We consider laboratory experiments that can detect stable, neutral strongly 

interacting massive particles (SIMPs). We explore the SIMP annihilation 



PsJ ' cross section from its minimum value (restricted by cosmological bounds) to 



the barn range, and vary the mass values from a GeV to a TeV. We calculate, 



^Sj . as a function of the SIMP-nucleon cross section, the minimum nucleon number 

''^ \ A for which there should be binding in a nucleus. We consider accelerator mass 

f~^ I spectrometry with a gold {A = 200) target, and compute the likely abundance 

On ' of anomalous gold nuclei if stable neutral SIMPs exist. We also consider 



the prospects and problems of detecting such particles at the Tevatron. We 
estimate optimistically that such detection might be possible for SIMPs with 
SIMP-nucleon cross sections larger than 0.1 millibarn and masses between 25 



Q ' and 50 GeV. 

>: 

•'-j ■ I. INTRODUCTION 

rS 
H , 
5t , The nature of the dark matter in the universe is a question of great interest for particle 

physics since the successful standard model has no candidate with the right properties to 

qualify as one. It is, therefore, hoped that determining the properties of the dark matter 

particle can provide important information about the nature of new physics beyond the 

standard model. One intriguing candidate that has been discussed in literature is a stable 

massive strongly interacting neutral particle (SIMP) |l||. Such particles arise in many gauge 

models |^, and their existence is of interest even if they do not constitute dark matter. 

While there already exist several severe constraints from cosmology and astrophysics on these 

particles, ||l|,|] they can still be viable if their strong interaction parameters are appropriately 

chosen. It is therefore of interest to explore all possible avenues to "tighten the noose" on 

them. In this paper we discuss two experiments of interest in this connection. 

The essential property of the SIMPs (henceforth denoted as X) that determines their 

possible experimental signature is the nature of their interaction with ordinary nucleons. It 

is not easy to guess the interactions of such unknown particles; however, we have some guide 

from nucleon-hyperon interactions. In addition, as detailed in Section II, we make use of 



an approximate relation between axN and <Jxx- Furthermore, cosmology provides a lower 
bound on ctj^y- In considering binding of X in nuclei, we will assume the XN potential to 
be attractive at low energies, as should be the case whenever XN scattering is a result of 
two particle exchange. 

In a recent publication |Q two of us (RNM and VLT) investigated constraints on SIMPs 
from searches for anomalous nuclei containing them. We considered two cases. First we 
built on the work of Dicus and Teplitz, ||^ assuming that the depth of binding potential 
VxN is of the order of a few MeV's, as is the case for the A particles. That work showed that 
the anomalous nucleus ^Be*x would be formed during big bang nucleosynthesis (BBN) if the 
X-N potential is attractive and large enough [i.e., lOMeV or more) for X-binding in light 
nuclei. Using a plausible model for the XX annihilation cross section, Ref. [Q estimated 
that the concentration of the X particles relative to the baryons must be bigger than ~ 10~^, 
which then leads to a concentration oi^Bex in the present universe which is several orders of 
magnitude higher than the present upper limits on anomalous Be obtained in experiments 
by Hemmick et al. 

Second, Reference |^ pointed out that, for a given attractive X-N potential, binding is 
more likely in the more massive high Z nuclei than in low Z nuclei. This is because, by the 
uncertainty principle, there is less kinetic energy for the SIMP in the high Z nucleus. We 
noted, however, the absence of data for anomalous nuclei with Z above that of sodium (A=23, 
Z=ll), and urged that experiments be undertaken in high Z nuclei. Since publication of 
Ref. 0, we have been informed P] that one group hopes to look for anomalous gold (A=200, 
Z=79) nuclei using accelerator mass spectrometry. 

The purpose of this paper is to build upon Ref [Q in two directions. First we compute 
X binding and abundance in gold as a function of Mx and axN- We choose gold in view 
of the experiment of Ref. [8]. Results for other heavy nuclei will be similar. Second, we 
make an initial estimate for detection of X-X pairs at the Fermilab Tevatron. The plan 
of the paper is as follows: in Section II we estimate the abundance of anomalous gold as a 
function of Mx and axN for the region of the {Mx, cfxn} plane for which there is binding. 
In Section III we consider Tevatron detection. Section IV gives our conclusions. 

II. SEARCHING IN GOLD 

We perform two calculations in this section. First, we compute X-A binding as a function 
of nucleon number A and Mx for a range of ctxn- Second, we compute the abundance of 
anomalous gold, relative to normal gold, as a function of Mx and axx- 

We consider first the nuclei in which the X-particle will bind. As in Ref. [^ and works 
cited therein, we take the X-A binding energy to be: 

where Vxn is the average X-N potential at rest, /x is the reduced mass of the X-A system, 
and R is the size of the nucleus in which X is bound. We have used a simple particle in a 
box model. We take i? ~ i?o A^f^ with Rq ~ 1.4 fermi, and A is the atomic number of the 
nucleus. The reduced mass of the X-A system is given by the usual expression: 
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The average X-N potential, Vxn, is defined implicitly via: 
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If we take the low energy NN scattering cross section a^^ to be on the order of 1 barn, 
use Vnn — 50MeV, and assume Bx <C Vx, we can write equation (|l]) as 
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which numerically reduces to 
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FIG. 1. The results of solving equation (|5|) for three choices for (Jxn- The dotted line is for 
a^-j^{min), the thin line is for c^xli ~ IM^X' ^'^d the thick line is for axN — <^nN- 

In Figure 1, we give the results of solving equation @ for three choices for axN- We 
determine the minimum value of A for which there is binding as a function of Mx- We 
assume an approximate factorization hypothesis of ref. M, 



P (^XX ^NN - Cr 



XN 



(6) 



where /5 ~ 1. 

To be definite, we shall consider three values for a 



XX- 



1. We consider cr^x ~ 3 x 10 ^^barn as a lower bound. This is the minimum cross section 
allowed which does not overdose the universe, (assuming Mx~ GeV). 

2. We use the estimate of ref. [|, a^x - 1/^1 = OAmb/Mj^{GeV) 



3. Finally we look at the estimate of ref. |^, axN — O.Ict^at ~ 2mb. This estimate takes 
X to be the particle responsible for ultra-high energy cosmic ray (UHECR) events 
0, the UHECRON. 



M 



For case (1), the minimum (J-^-j^, we see in Figure 1 that gold will only capture X's for 

1/M|^, we estimate 



x^80, (assuming the interaction is attractive). For case (2), a-^-j^ 



there is always binding in relatively light nuclei, In particular, much of the Mx range (15~150 
GeV, roughly) would have binding in sodium, and hence is likely ruled out by current 
experiments. For case (3), the UHECRON, there is binding in gold for Mx above a few 
GeV. 

We turn now to the less straightforward problem of estimating the abundance of anoma- 
lous gold. Our first requirement is a capture scenario. Two possibilities present themselves. 
First, the gold on Earth was likely made by a supernova in the giant molecular cloud (GMC) 
in which the sun was bound. GMC's live for tens of million of years, so one could assume 
capture in the GMC before formation with an exposure time on the order of 10^ years. 
There are complications, however. The X needs to be slowed from the galactic virial veloc- 
ity down to thermal velocities, and then captured. In addition, both the capture reaction, 
X + A —* ^ A + 7, and the dissociation reaction, 7 -|- ^ A -^ X + A from 7's generated by 
hot stars and supernova in the GMC, must be considered. There appear to be too many 
unknown parameters for a reliable estimate. The second, more conservative and more easily 
estimated scenario is to assume SIMP capture after formation of the Earth. 

The key to the second scenario is that gold nuclei be close enough to the surface for 
galactic halo SIMPs to reach them. Geologists say, with some confidence, [jll| that there are 



selected gold deposits that should have had long exposure to any cosmic SIMPs. These are 
deposits found at the surface of the earth, in particular, in the gravel of streams ( "placer 
gold") in regions that are sufficiently inactive geologically that one can be reasonably con- 
fident of an exposure for over ~ 10^ years. The Sierra Nevada Mountains are such a region. 
Placer gold from the California gold rush would have a lesser exposure time, on the or- 
der of ~ 10^ years. For such placer gold, the results discussed below on the abundance of 
anomalous gold should be decreased by a factor of 10^. 

Given an exposure time of ~ 10^ years, we must then calculate how deep into the Earth 
the X particle penetrates, and how many nuclei are in the region penetrated. We make the 
conservative approximation that the gold is uniformly mixed with other elements so X can 
be captured by any A > Amin- We assume all X's are captured since scattering slows them 
to thermal velocities. 
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Element 


Cm in {figm/gm) 


< 20 


20 


Ar 


5x 10^ 


10 < A < 26 


26 


Al 


8.3 X 10^ 


26 < A < 137 


137 


Ba 


425 


137 < A < 200 


207 


Ph 


12.5 



TABLE I. Abundances of the most abundant nuclei in a given A-interval. 



We now compute the fraction / of gold that has captured an X: 

f{Mx,crxN)^nxvt^ (7) 

Here nx is the galactic halo abundance at the solar distance from the galactic center, v is 
the galactic virial velocity (300 km), t is the exposure time, /au is the fraction of the stopped 
X's captured in gold, and Nau is the number of Au nuclei (per cm^) in the stopping length 
Ax of X. 

The fraction Jau is taken to be: 



Cau J Am . s 

where the C's are concentrations by weight in the Earth's crust, and we expect the capture 
cross to be proportional to A^^^. The index M designates the most abundant element with 
a nucleus heavier than A^in- We choose Am from Table I. 

Nau in equation (|^, the number of Au nuclei (per cm^) in the stopping length Ax of X, 
is given by 

Nau = Ax ^ ^ (9) 

where Me is the mass of the earth, and nip ~ 1 GeV is the mass of an ordinary nucleon. 
For Ax we take 



Vx = S T^ A'/' axN 
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where A ~ 20 is the A value for an "average nucleus," and the last factor (= AEx/Ex) 
assumes that XA scattering in the CM. is isotropic. Note, as a sanity check, that the Cau 
factor in equation(|) cancels with that in /au, equation(H), as one would expect. 

Finally, to get nx in equation(H) we compute the primordial X freeze-out abundance, 
nx, precisely as specified in Kolb and Turner |T^, using (in their notation) n=l and g=l, 



and making use of the factorization approximation relating Cxx ^'^ (^nx- We then assume 
that X's are concentrated in the galaxy to the same extent as baryons, but are distributed 
in an isothermal spherical halo of 70kpc radius. This gives nx/ns = 3.7 x lO^^fix/nB with, 
the tilde denoting the cosmic average. 

The result of the calculation for logio(l//) is shown with quiet drama in Table 2. axN 
ranges from 5 x 10^^6 (corresponding to a-^j^ = amin = 3 x lO^^^cm^) to 1 barn, while Mx 
ranges from 1 GeV to 1 TeV. In the first column we see that the X does not bind in gold until 
Mx e [80, 130] GeV, which is consistent with Figure 1. As axN increases, the abundance 
falls because [0 Yoo oc l/{Mxcrxx) &nd aj^-j^ ~ (^nx- This fact also plays a major role in 
the decrease of / as Mx increases. Note from the table that the boundary between binding 
and no-binding can be approximated by the implicit relation MxaxN ~ bmbGeV'^. 

The result of Table 2 is optimistic in that previous experiments have reached abundance 
fractions approaching / ~ 10^^°. Therefore, it is likely to be possible to explore the entirety 



of the parameter space of Table 2 in which there is binding of the X. For the UHECR case 
of Ref. H (ax7v~0.1or^jv, M^^SO GeV), it is only necessary to go to IQ-^^ not IQ-^o. This 
is because Table 2 shows that, even for the largest XN cross sections, if the value of Mx is 
below 50 GeV, logio(/) > -15.7. 

While Table 2 shows the fraction / is a relatively smoothly varying function of a^x and 
Mx over most of the range, closer examination shows some regions of large variation. For 
example, the first column shows a change by 2 orders of magnitude between Mx=350 and 
570 GeV. In Figure 1 this corresponds to Mx — 530 where A^i^ goes from above 137 to 
below 137. This corresponds to Am in equation(^) falling from 207 (Pb) to 137 (Ba) which 
shifts Cm in equation(|]) by a factor of 35. To see this effect more clearly, we plot in Figure 2 
abundance fraction / as a function of Mx for the 3 cross section choices of Figure 1. The 
curve for a-f^-j^{min) has a steep rise just past Mx = 100 where Am in Figure 2 first falls 
below 200 (so that there is binding in gold). The discontinuous behavior at Mx = 530 GeV 
is seen clearly. The curve for axN — (^-kN is monotonic (after its initial rise when binding 
occurs) corresponding to the monotonic behavior of Au in Figure 1 for the corresponding 
curve. The curve for c^xjc — ^/^x has discontinuities of both signs corresponding to the fall 
of this same curve in Figure 1 to below 26 and 20, and subsequent rise above these points. 
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TABLE II. Mx (vertical) is in units of GeV, and ctxn (horizontal) is in units of mb. Table 
entries are log]^o(l//)i ^^^ ^^^ ~ indicates those cases for which X does not bind at all. 
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FIG. 2. The abundance fraction / as a function of Mx for the 3 cross section choices. The 
dotted hue is for a-^j^{min), the thin hne is for cr-^j^ ~ 1/Mj^, and the thick line is for axN — cTwN- 

It is also of interest to make a connection with the experiments searching for weakly 
interacting massive particles (WIMPs). For crxAf~4 x 10~^6, which corresponds to a^^x ^ 
factor of 50 above o"j^^(min), the contribution of X to the local galactic dark matter density 
(4 X 10~'^^g/cm^) is less than 2 percent. 

It is important to note that we only assume that the mass density of SIMPs is less than 
or equal to the dark matter constraint; we do not assume, as did Starkman et al, |l[], that 
it saturates the constraint. 

In summary, we see from the above discussion that: i) we have binding in gold for 
MxaxN^^mbGeV'^; ii) a SIMP of mass up to a TeV satisfying this condition, and having 
an attractive interaction with nucleons, should form anomalous gold nuclei of sufficient 
abundance to be seen by accelerator mass spectrometer experiments sensitive to one part in 
10^°; and iii) UHECRs could be detected in an experiment with a sensitivity of one part in 
10^^ for the mass range from a few GeV to their 50 GeV maximum, (assuming a target of 
gold with 10^ years' exposure). 



III. SEARCHING AT THE TEVATRON 



If neutral, stable SIMPs were actually to exist, it might be possible to produce and detect 
them at the Tevatron. Although such detection might be difficult, further study is required 
before one can decide whether it is impossible. In this section, we make an optimistic 
estimate of the most promising ranges of Mx and axN for observing these SIMPs. While 
the region accessible to the Tevatron is quite restricted, it includes a portion of a particularly 
interesting region - the region relevant to the explanation of the Ultra High Energy Cosmic 
Ray (UHECRs) events proposed by Farrar, Kolb and coworkers. |^ 

We consider XX production in one of the Tevatron detectors. Assuming that the SIMP 
has colored constituents {e.g., gluon, gluino, ...), its pair production will be accompanied by 
soft hadrons. A neutral, stable SIMP particle will not generate a signal in the central tracker, 
and is unlikely to interact in the electromagnetic calorimeter of the Tevatron detectors. It 
may, however, interact with material in the denser hadronic calorimeter. Such an interaction 



will have kinematical characteristics that can be used to distinguish it on a statistical basis 
from neutron and K^ interactions. The distinction is based on the assumption that at 
the Tevatron energies, the mass of a SIMP represents a substantial fraction of its total 
energy. We then need to ask how much energy the SIMP is likely to deposit in the hadron 
calorimeter, whether we can distinguish a SIMP shower from that from a neutron or K^, 
and what overall event rate might be expected. 
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FIG. 3. Contours for energy loss as a function of {Mx-,Ex}- The contours displayed are in 
steps of 10 GeV. 

With regard to energy deposition in the hadron calorimeter, we make a qualitative es- 
timate. Assume that, in interacting with a nucleon in one of the steel plates, the X loses 
half of its CM. kinetic energy into particle production. We can then transform back to 
the lab system and estimate the energy lost by X as a function of Ex and Mx- Figure 3 
shows the results. It has contours for energy losses (into hadron showers) of 10, 20,... 100 
GeV. We see that 400 GeV partons could, triggering on multiple ~10 GeV showers, permit 
discovery of SIMPs with masses up to 60 GeV. Similarly 600 GeV partons could push the 
SIMP discovery limit up to 100 GeV. At 300 GeV, we can get up to 50 GeV masses, which is 
the upper limit, in the analysis of Farrar et al 0, of a SIMP that would explain the UHECR 
events. 

If axN is too small, there will be no showers at all. In Im of steel, we get about 
5 interactions for a cross section of a few m6, but only 0.001 for our minimum axN (~ 
5 X lO^^^cr/i^). We might detect one jet from the SIMP pair at a sufficient rate if we had a 
cross section on the order of at least a tenth of a millibarn, and both jets for a cross section 
over a millibarn. The estimate of reference states that a cross section of more than a 
millibarn is needed if the SIMP is to reproduce the UHECR events. 

We estimate the production rate of SIMPs by scaling the production rate of jets by the 
ratio of XN to the Meson-N cross sections. We observe that the total cross section for vrp 
and Kp scattering is approximately ~ 30m6 for ^/s of a few hundred GeV, with a mild 
(logarithmic) variation with ^/s. With this perspective, we can now turn from the t-channel 
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Xp —>■ Xp process to the s-channel pp — > XX to estimate a scaled cross section. In the low 
energy region, X-production will be suppressed relative to jet production by the X mass. 
However, for Ex larger than a few times Mx, there is no more phase space suppression, and 
we therefore expect scaling to work well (for Mx = 50 GeV) in the region Ex^'200GeV. 
The cross section for producing jets with ET^200GeV is ~ lOOpb. This gives ~ 3pb for any 
one quark (to be compared with a t-quark production cross section of ~ 5pb) , corresponding 
to ~300 events in Run I and ~6000 events in Run II. Again, reference estimates that 
the UHECR would need a cross section of about a tenth the meson nucleon cross section. 
Scaling the production cross section by that amount gives a good number of events in Run II. 
Turning to the question of whether we would recognize a SIMP shower in a hadron 
calorimeter if we saw one, the best available discriminant would appear to be the opening 
angle of the shower. A pion moving transverse in the SIMP-nucleon CM. system will have 
a laboratory angle of 



-f Ex + nip 

Comparing a 500 GeV SIMP with a 500 GeV neutron, we see that the SIMP shower will be 
~ 40 percent wider if Mx^25GeV. 

In summary, an optimistic scenario finds the Tevatron discovery potential in the range 
CTxAT^O.lmfe and 25GeV^Mx^50GeV. Consequently, it might be possible to resolve the 
UHECR puzzle at Fermilab. 



IV. CONCLUSIONS 

We emphasize here that, as discussed in Ref. [^, searches for SIMPs in anomalous nuclei 
are much better carried out in high Z nuclei. In such nuclei, the SIMP kinetic energy is 
minimized so that the chances of binding are maximized. Table 2 shows that a proposed 
accelerator mass spectrometry experiment in gold is capable of discovering or eliminating 
SIMPs over ~80 percent of the relevant portion of the {(Txn, Mx) plane so long as the low 
energy X-N potential is attractive, (the likely event). 

Looking for X at the Tevatron is difficult. One needs events in which there are no high 
energy particles at the vertex or in the electromagnetic calorimeter, but there are showers in 
the hadron calorimeter. The range of Mx values for which the hadron calorimeter showers 
are sufficiently energetic goes up to 50 or 100 GeV for the Mx mass. The range for which 
the shower opening angles are sufficiently wide to distinguish them from neutron or Kaon 
showers begins around 25 GeV for the Mx mass. Although this region is limited, it includes 
at least half the region of interest for a possible UHECR explanation. |^ 
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